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ABSTRACT 

We conduct a search for radio stars by combining radio and optical data from the Faint Images of 
the Radio Sky at Twenty cm survey (FIRST) and the Sloan Digital Sky Survey (SDSS). The faint limit 
of SDSS makes possible a homogeneous search for radio emission from stars of low optical luminosity. 
We select a sample of 112 candidate radio stars in the magnitude range 15 < z < 19.1 and with radio 
flux 5*20 > 1.25mJy, from about 7000 deg^ of sky. The selection criteria are positional coincidence within 
1", radio and optical point source morphology, and an SDSS spectrum classified as stellar. The sample 
contamination is estimated by random matching to be 108 ±13, suggesting that at most a small fraction of 
the selected candidates are genuine radio stars. Therefore, we rule out a very rare population of extremely 
radio-loud stars: no more than 1.2 of every million stars in the magnitude range 15 < i < 19.1 stars has 
radio flux 820 > 1.25 mJy. We investigate the optical and radio colors of the sample to find candidates 
that show the largest likelihood of being real radio stars. The significant outliers from the stellar locus, as 
well as the magnetically active stars, are the best candidates for follow-up radio observations. We conclude 
that, while the present wide-area radio surveys are not sensitive enough to provide homogeneous samples 
of the extremely rare radio stars, upcoming surveys which exploit the great sensitivity of current and 
planned telescopes do have sufficient sensitivity and will allow the properties of this class of object to be 
investigated in detail. 

Subject headings: stars: statistics — radio continuum: stars — surveys 
1. INTRODUCTION 
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processes that lead to these types of emission — plasma 
heating and particle acceleration in stellar coronae — 
are seen in our own Sun, but the relevant energies for 
active stars are much larger. Radio emission at the 
relative level of that emitted by the Sun remains unde- 
tected from even the closest solar-type main sequence 
stars to the present day0 The quiescent, slowl y vary- 
ing ra dio emission seen in many active stars (e.g. Giidell 



2OO2L and references therein) has no solar counterpart. 



With the great increase in the sensitivity of ra- 
dio surveys in the last several decades, along with 
the more accurate source positions allowed by radio 
interferometry, both thermal and non-thermal radio 
emission have now been detected from hundreds of 



stars of many differ ent tvpes faiellming fc Gibson 
1986[ IWendked [l987t lAltenhoff et al.l ll994l: IWendkei 
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1993 ). rapidly - rotati ng main-sequence stars (jLim fc White! 
1995t Bergei 2002 ). X-ray bright ma in sequence 



stars (Giidel et all 1995 ). magnetic stars dPrake et al 
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1986 : lDrake et al.lll987bl : lKnapp et al. 1995; Reid fc Menten 
I997I) OB stars wi th winds dBieging et all Il989[ 
Drakdll99d [P hillips fcTitiijll990l), Wol f-Rayet stars 



..Chapm an et all 119991). dMe flare stars (jWhite et al, 
19891 : lOsten et al.ll2006f ). and various classes of inter- 



acting binaries and cataclysmic variables. The radio 
radiation from these stars is very faint, at the few mJy 
level. 

The first large unbiased study of radio stars (~ 
5000 deg^ of sky at high galac tic latitude) was per- 
formed by iHelfand et all |l999l . hereafter H99), who 
compared the Faint Images of the Radio Sky at Twenty 
cm survey (FIRST; iBecker et al.lll995^ to several cat- 
alogs of bright stars with high astrometric precision: 
the Hipparc os catalog (Ferryman et al. 1997), the Ty- 
cho catalog (iHoeg et al.lll997t ). the Guide Star Cata- 
log ([Lasker et al.lll990l ). and stars within 25 pc of the 
Sun. They emphasize the need for accurate positions: 
the rarity of radio stars to the FIRST flux limit (~ 
1 mJy), combined with the high density of faint ex- 
tragalactic radio sources, ensures random matches be- 
tween stellar and radio sources in sufficient numbers to 
confuse the cataloging of true radio stars, unless both 
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radio and optical positions are known to better than 
1". H99 identified 26 radio stars in their study, about 
one per 190 deg^, and showed that the fraction of stars 
with radio emission above the FIRST limit declines 
steeply with optical magnitude to my ^15. The frac- 
tion of radio star s at fa inter magnitudes is unknown. 
Kimball fc Ivezid (|2008l hereafter KI08) searched for 
radio stars in a combined radio — optical catalog with 
observations from FIRST and from the optical Sloan 
Digital Sky Survey (SDSS). Quasars are the most com- 
mon radio source above flux densities of a few mJy; 
therefore a sample of optical point sources with radio 
emission is likely to be srongly dominated by quasars. 
KI08 approached this problem by applying a conserva- 
tive photometric color cut, using the fact that quasars 
and stars lie in diff erent locations in SD SS optical 
color-color diagrams (iRichards et al. I I2OOII) . The sam- 
ple was limited to sources sufficiently bright to be in- 
cluded in the SDSS quasar spectros copic target selec- 
tion {i < 19.1. ISchneider et alllioOTh . Only 20% of the 
photometrically-selected candidate radio stars actually 
showed stellar spectra, while the rest were quasars with 
stellar-like colors. KI08 concluded that simple color 
critera are not sufficient to select a clean sample of 
radio stars, and that spectroscopic observations are 
necessary to distinguish between quasars and stars. 

In this paper we continue the search for radio stars 
in the SDSS using a sample with spectroscopic identifi- 
cations. We present 112 candidates selected by match- 
ing FIRST detections and SDSS point sources within 
1 arcsec. The sample comprises sources brighter than 
r = 20.5 in the optical and 1.25 mJy at 20 cm, with 
SDSS spectra classified as stellar both by the auto- 
mated reduction pipelines and visually. The SDSS 
spectroscopic targeting implicitly imposes soft mag- 
nitude limits of 15 < i < 19.1. In this magnitude 
range, approximately 1% of SDSS stars have spectro- 
scopic data. However, all objects in this range which 
are close to a FIRST source are targeted for spectra 
(jStoughton et al.l 120021 ). so the completeness of SDSS 
radio — optical sources is well understood. 

We are searching in a different region of the radio — 
optical parameter space from II99, who also matched 
to FIRST but used optical catalogs brighter than 
TTiy = 15. The SDSS has a saturation limit of = 15. 
By extending to several optical magnitudes fainter 
than H99, we are therefore searching for stars with 
a much brighter radio-to-optical flux ratio. However, 
a FIRST — SDSS matching has the potential to re- 
veal a radio-bright population too rare to appear in 
the smaller H99 study, as the fainter SDSS includes 
a much larger volume in the Galaxy. For example: 
M dwarfs, many of which are kno wn to be magneti- 
cally active (e.g., West et allbOOSl) . are found in sig- 
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nificant numbers only at faint magnitudes. Thus, a 
search for radio-emitting M dwarfs should be carried 
out in a deep optical survey with large sky coverage. 
Advantages to using the SDSS are its high complete- 
ness, precise magnitudes, accurate astrometry, much 
deeper optical data than previous stellar catalogs, and 
almost 300,000 stellar spectra. 

The remainder of the paper is laid out as follows. 
In 31] we describe the contributing surveys. In iJ3] we 
outline the selection of the sample of radio stars and 
place an upper limit on the fraction of radio stars. In 
§4] and 50 we discuss optical and radio properties of 
the sample, respectively, and in ^ we conclude and 
summarize our results. 

2. OPTICAL AND RADIO DATA 
2.1. Optical Catalog: SDSS 

We have drawn our sample from the photomet- 
ric coverage of the sixth data release (D R6) of the 
Sloan Digital Sky Surve#l (SPSS: see lYork et al " 



(jl998r ). The ast rometry of the p hotometric survey is 



200: 



2000l:IStoughton et al.ir2002l : lAdelman-McCarthv et al 



2008L and references therein). DR6 covers roughly 
9, 600 deg^ and contains photometric observations for 
287 million unique objects, as well as spectra for 
more than 1 million sources (in a smaller sky area 
of 6860 deg^). SDSS entered routine operations in 
2000; DR6 observations were completed in June 2006. 
Because SDSS spectroscopy is performed after pho- 
tometric observations, some stars in the sample have 
spectra that were not available until the Seventh Data 
Release (DR7). We have included some stars with 
DR7 spectrsEI; we point out explicitly where this de- 
tail affects our estimates of sample contamination in 
§3.61 The DR6 sample includes about 287,000 spectra 
classified as stars. 

The SDSS photometric survey measures flux den- 
sities nearly simultaneously in five wavelength bands 



4686, 6165 


7481, and 893lA (Fukueita et all 


1996 


Gunn et al. 


1998; HoEK et al.ll2001;lTucker et al. 


2006 


Smith et al. 


2OO2I; Ivezic et al.l 20041). Morohologv in- 



formation allows reliabl e star - 

2OO2I : 



21.5 (jLupton et al 



galaxy separation t o 



Scranton et"all l2002h . 



Sources are classified as resolved or unresolved using 
a measure of light concentration that determines how 



well t he flux resembles a point source (jStoughton et al 
20021 ). Magnitudes were corrected fo r Galactic ex- 



tinction according to the dust maps of iSchlegel et al 



^Tho survey website is located at http://www.sdss.org. 
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5 Dec 2007. 



good to < 0.1" (jPier et al. 

A subset of photometric sources was chosen for 
spectroscopy according to the SDSS targeting pipehne; 
stars could be serendipitously selected by any of the 
various targeting algorithms. The quasar targeting al- 
gorithm selects all 15 < i < 19.1 sources within 2" 
of a FIRST catalog object, and some sources as faint 
as i — 20.5 (depending on the availability of spectral 
flbers). About 30% of quasar targets turn out to be 
stars or galaxies ( Schneider et al.ir2007t ) . Another algo- 
rithm targets interesting stellar classes by selecting for 
their distinctive photometric colors; these include blue 
horizontal-branch stars, carbon stars, subdwarfs, cat- 
aclysmic variables, brown and red dwarfs, and white 
dwarfs. Because stars in these catagories are selected 
randomly to fill excess spectral fibers, the complete- 
ness of the samples of stars with spectroscopic obser- 
vations is, with a small number of exceptions, not well- 
defined. However, the fact that all objects within 2" of 
a FIRST source are targeted implies that the spectro- 
scopic sample is complete to i < 19.1 with respect to 
radio — optical sources brighter than the FIRST limit, 
except in the case of fiber collisions (see below). 

SDSS spectra are obtained using 3" fibers in a fiber- 
fed spectrograph (|NewmBiL£LalJl200j). Spectra cover 
the wavelength range from about 3500 to 9500 A with 
a resolution of R^^ISOO. They are automatically ex- 
tracted and calibrated by the Spectroscopic Pipeline: 
wavelength calibrations are calculated from arc and 
night sky lines and flux calibrations from observations 
of standard F stars. Some fibers are used to observe 
blank sky in order to correct targeted objects for the 
sky spectrum. Not every object targeted for spec- 
troscopy obtains a spectrum, primarily because of fiber 
collisions: due to the physical size of the fibers, no two 
fibers can be placed closer than 55" on a spectroscopic 
plate. As a result, spectra are only taken for about 
92.5% of spectroscopic targets, although the complete- 
ness effects are well understood (jBlanton et al..,2003 ). 

2.2. Radio Catalogs 

FIRST provides 20 cm radio fluxes, which were used 
to select candidate radio stars. The Westerbork North- 
ern Sky Survey and the Green Bank 6 cm survey pro- 
vide fluxes at 92 cm and 6 cm, respectively, which 
allow us to determine a radio spectral shape. 

2.2.1. FIRST 

The Faint Imag es of the Radio S ky at Twenty 
cm survey (FIRST; iBecker et aDll995l ) used the Very 
Large Array to observe the sky at 20 cm (1.4 GHz) with 
a beam size of 5"4 and an rms sensitivity of about 
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0.15 mjy beam ^. Designed to cover the same region 
of sky as the SDSS, FIRST observed 9, 000 deg^ at the 
North Galactic Cap and a smaUer ^ 2.5° wide strip 
along the Celestial Equator from 1994 to 2002. The 
survey contains over 800,000 unique sources, with posi- 
tions determined to < 1"; its source density is roughly 
97deg~^. It is 95% complete to 2mJy and 80% com- 
plete to the survey limit of 1 mJy. The integrated 
flux density, S20 , is calculated using a two-dimensional 
Gaussian fit to each source image. 

2.2.2. WENSS 

The Westerbork N orthern Sky Survey (WENSS; 



Rengelink et al.l 119971 ) is a 92 cm radio survey that 
was conducted in the mid-1990s. It observed the sky 
north of 6 — 29° to a limiting flux density of 18 mJy, 
with a beam size of 54" and a positional uncertainty 
of 1.5-5". 

2.2.3. GB6 

T he GB6 survey at 4.8 5 GHz (Green Bank 6 cm sur- 
vey; Gregorv et al. 19961 ) was executed with the origi- 
nal 91m Green Bank telescope in 1986 November and 
1987 October. Data from both epochs were assembled 
into a survey covering the 0° < 5 < 75° sky down to 
a limiting flux of 18mJy, with 3.5' resolution and a 
positional uncertainty of 10-50". 

3. SELECTING A SAMPLE OF CANDI- 
DATE RADIO STARS FROM FIRST 
AND SDSS 

This section outlines the selection criteria that yield 
the candidate radio stars sample. We estimate the 
amount of contamination with random matching, and 
use the result to place an upper limit on the fraction 
of radio stars in the magnitude range of the sample. 
The area of overlap of the two surveys is about 9500 
deg2. 

3.1. Magnitude and flux limits 

We applied magnitude/flux limits as a way to con- 
trol source quality. For SDSS sources, we required 
r < 20.5 to ensure reliable determination of optical 
morphology, as well as to select sources with a high 
enough signal-to- noise ratio (S /N) for spectral typing. 
Owing to the magnitude limits of the SDSS quasar 
target selection algorithm ( ij2.ip . most of the sources 
in the final sample have i < 19.1, although it is not 
an explicit requirement. There are roughly 4000-5000 
SDSS sources with r < 20.5 per square degree on the 
sky, depending on Galactic latitude. 



For the radio sources, we adopted a limit of 6*20 > 
1.25 mJy (equivalent to AB magnitude 16.2). This is 
slightly brighter than the 1 mJy depth of the FIRST 
catalog; however, in a visual examination of faint 
FIRST images, many sources fainter than this ap- 
peared to be possibly spurious detections. Above this 
flux limit, there are about 82 FIRST sources per square 
degree. 

3.2. Positional matching of point sources 

We restricted the sample to optical point sources 
using SDSS automated star — galaxy separation ( §2.ip . 
Figure [1] shows the distribution of FIRST— SDSS dis- 
tance for radio — optical (point source) matches. The 
expected contamination by random matches (evalu- 
ated by off-setting the FIRST positions by 1° in right 
ascension) is also shown. The inset plot shows the esti- 
mated completeness (solid line; percentage of physical 
sources recovered) and efficiency (dotted line; percent- 
age of all matches which are physical) as a function 
of matching radius. The precise choice of matching 
radius is a tradeoff between sample completeness and 
sample contamination. Since we are looking for rare 
objects, we opted more on the side of decreased con- 
tamination, and chose to use a 1" matching radius. As 
shown in Figure [H the matching radius of 1" results 
in 96% completeness and 98% efficiency. There are 
2000-3000 SDSS point sources per square degree with 
r < 20.5. Correlating with FIRST positions within 1" 
resulted in 14, 000 matches. 

We did not consider proper motions when match- 
ing the two catalogs, but emphasize that this decision 
should not significantly affect the matching results. All 
of the observations were performed since 1994. Within 
the relevant magnitude range, nearly 99% of SDSS 
stars have an apparent motion of less than 0.1 arc- 
sec yr~l, with a median value < 7 mas yr~^. Proper 
motionqj for this sample are much smaller than for the 
sample of H99 owing to the much fainter flux limit and 
correspondingly larger source distances. 

3.3. Spectral typing 

KI08 showed that a sample of potential radio stars 
selected by their SDSS colors is strongly contaminated 
by quasars. There are enough quasars with stellar col- 
ors (i.e., within a few hundredths of a magnitude from 
the main stellar locus in the multi-dimensional SDSS 
color space) that spectroscopic identiflcations are nec- 
essary in order to cull them from the sample. We 
therefore required an SDSS spectroscopic observation 



^Proper motions were determined from t he Seventh Dat a Release 
of SDSS following the procedure of ,Munn et al.l ||2004D . 
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for each matched source, and hmited the sample to 
those with stellar-classified spectra. Of the ~ 14, 000 
selected matches, 6413 have spectra, and 292 of those 
were classified as stellar by the SDSS spectral process- 
ing pipeline. 

For a more robust determination of spectral type, 
we executed visual classification of each spectrum. Of 
the 292 stellar-classified sources, 6 matched to known 



quasars (ISchneider et al.l 120071) a nd 45 to known BL 



Lac objects ( Plotkin et al. 20081 ). We removed one 



object whose spectrum showed broad emission lines 
signifying the super-position of a star with a quasar. 
Of course these false matches are not indicative of an 
18% failure rate for the spectral pipeline, which is op- 
timized to detect galaxies, stars, and quasars rather 
than rare sources such as BL Lac objects. Because the 
selection criteria for this study are biased toward (ra- 
dio) quasars, they result in a disproportionately large 
number of BL Lac objects, which make up a very tiny 
percentage of the entire SDSS spectroscopic database. 

Spectral types were individually assigned to the 
remaining 240 sources using a custom IDL pack- 
age dubbed "the Hammer" |j The full algorithm 
us ed by the H ammer is described in Appendix A 
of ICovev et al.l ([2007). In short, the Hammer auto- 
matically types input spectra by measuring a suite of 
spectral indices and performing a least-squares min- 
imization of the residuals between the indices of the 
target and those measured from spectral type stan- 
dards. It then allows a confirmation or correction of 
spectral type according to a visual comparison of the 
input spectrum with spectral templates. Although 
spectral types are available from the SDSS database, 
visual confirmation by stellar scientists (authors GRK, 
AAW, JJB) leads to more robust classifications. 

With visual classification, the sample was reduced 
to 194 sources with reliable stellar spectra. The re- 
jected spectra were too noisy for reliable typing, or 
were indicative of BL Lac objec ts. The latter are not 
included in the BL Lac sample o f lPlotkin et alllEoOH) : 
that sample was drawn from SDSS Data Release Five 
(DR5), and thus covers a smaller sky area than the 
DR6 sample used in this paper. 

3.4. Visual examination of radio morphology 

We looked at FIRST images of the remaining sam- 
ple in order to determine the radio morphology of each 
source. We anticipated that the majority would be 
point sources, because the resolution of FIRST (~ 5") 
is not sufficient to resolve stellar emission. Resolved 
or multiple-component emission, however, would be 



strongly indicative of a non-stellar source such as an 
active galactic nucleus (AGN) with radio jets. 

We examined FIRST postage stamps (2' x 2'), 
and classified each as "compact" (point source emis- 
sion), "resolved" (resolved single-component emis- 
sion) , "complex" (multiple-component or knotty emis- 
sion), or "spurious" (e.g., an artifact introduced by in- 
terferometric errors). An example from each category 
is shown in Figure [51 Resolved or complex radio mor- 
phology is typically associated with an extra-Galactic 
object such as a radio galaxy or quasarQ Out of the 
194 visually-classified images we found six spurious 
sources, which were removed from the sample. Of 
the remaining 188 sources, 60 are complex, 16 are 
resolved, and 112 are compact H The high fraction 
of sources with complex radio emission demonstrates 
that many radio quasars survived the previous selec- 
tion criteria. Because these objects clearly show stel- 
lar spectra, we interpret them as optically-faint radio 
quasars in chance alignment with bright foreground 
stars. We rejected the complex sources as obviously 
extra-Galactic. However, not all quasars have complex 
radio emission from detectable lobes, but many are in- 
stead point sources (KI08): this is typically thought 
to be the result of Doppler beaming of a jet aligned 
along the line of sight. Thus, despite the thorough 
visual classifications, the final set of compact radio 
sources with stellar spectra may remain contaminated 
by chance star — quasar superpositions. We discuss 
sample contamination in more detail in 



3.5. Final sample of potential radio stars 

We retain the sample of 112 compact sources for 
the remaining analysis of this paper. This set of radio 
star candidates is presented in Table [TJ which lists po- 
sitions, stellar type, radio fiuxes, optical magnitudes, 
and radio morphology classification. The 16 resolved 
sources and 60 complex sources are listed in Tables [2] 
and [HI respectively. The tables also list the distance 
to each sta r, calculated using th e photometric parallax 
relation of llvezic et al. ( 2008a ): this relation is valid 
only for stars o n the main sequence. As discussed by 



Finlator et all (|200fl[ ). nearly aU SDSS stars (~ 99%) 
are expected to lie on the main sequence. However, 



^http://w'ww.cfa. harvard, edu/^kcovey/thehammer 



^Galactic objects such as HII regions or supernova remnants may 
show extended emission; however we do not expect contamina- 
tion from these sources because our sample is hmited to 30° or 
more above the Galactic plane. 

*For one source that was initially identified as "complex", the 
optical image shows a nearby galaxy that appears to be the 
source of the second point of radio emission. We concluded that 
the two radio components are physically unrelated, and moved 
that object into the "compact" category, loading to the totals 
given in the text. 
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it is possible that by selecting radio-emitting stars we 
have biased our sample toward giant stars in the halo 
(nearby giants are brighter than the to = 15 SDSS 
saturation Hmit). We note that approximately 40% 
of the candidate radio stars have \og{g) measured by 
the SDSS software pipeline; all have log(g) > 3, in- 
dicating that they are main sequence stars. There is 
no evidence that the sample has a different distance 
distribution than other SDSS stars with spectra. 

Figure [3] is a mosaic of four of the candidate radio 
sources; it shows the optical image, the radio image, 
and the optical spectrum for each one. The brightest 
optical source and the brightest radio source in the 
sample are included in the mosaic. 

3.6. Estimating Sample Contamination 

To estimate the contamination originating from 
chance radio — optical alignments, we created a set of 
twelve random samples for comparison. To create the 
random samples, we off-set the right ascension or dec- 
lination in the FIRST catalog (by -1, -0.5, -0.1, 0.1, 
0.5, or 1 degree), then apphed identical selection cri- 
teria (where possible) as above. The contamination 
estimate is an upper limit because it is not possible to 
apply exactly the same selection criteria to the ran- 
dom samples: the SDSS spectral target selection for 
quasars ( §2.ip ensures that a large fraction of the real 
FIRST — SDSS matches have spectral data, whereas 
very few sources in the random samples have SDSS 
spectra. 

Having performed the matching several times, we 
can determine the variance of the random sampling, [f] 
Random matching of all FIRST and photometric SDSS 
sources resulted in 3242 ± 65 matches. Applying the r 
magnitude and radio flux limits reduced the samples 
to 759 ± 30 matches. Selecting on optical point source 
morphology further reduced the samples to 422 ± 20. 

The next selection step for the real sample was to 
eliminate those sources without SDSS spectra. The 
equivalent step for the random samples is to eliminate 
those which could not have qualified for SDSS spectral 
targeting. The quasar targeting algorithm ( §2.ip de- 
pends on proximity to a radio source; by artificially 
shifting the FIRST positions to create the random 
samples, we created fake optical — radio sources which 
would pass the selection criteria. We rejected sources 
outside of the DR6 spectroscopic coverage, which is 
smaller than the DR6 photometric coverage. Keep- 



^ Varying the FIRST positions results in a slight decrease of the 
areal overlap between FIRST and SDSS. However, given that so 
few matches result from the random sampling, the small change 
in matching area has no significant effect on the contamination 
estimate. 
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FIRST-SDSS distance ["] 

Fig. 1.— The distribution of FIRST— SDSS distances 
for FIRST sources matched to SDSS r < 21 point 
sources (diamonds). The level of background contam- 
ination, (estimated by offsetting positions by 1°) is 
also shown (squares). The inset plot shows the com- 
pleteness (solid line) and efficiency (dotted line) as a 
function of matching radius, estimated by fitting the 
distance distribution with a simple Gaussian + line 
model. 





• 







Fig. 2.— Example FIRST images (2' by 2') classified 
by radio morphology as spurious (upper left), compact 
(upper right), resolved (lower left), and complex (lower 
right). The position of the optical source is at the 
center of each image. 
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Fig. 3.— Data for four of the FIRST— SDSS star 
matches. The left column shows the optical image 
(2'x2'), a composite of the 5, r, and ? band images; the 
middle column shows the FIRST image (2' x 2') with 
linear scaling; the right column shows the SDSS spec- 
trum (Fx). The top row shows the brightest optical 
source, an Ml dwarf with i = 14.9 and S20 = 4.3 mJy. 
The second row shows a G3 dwarf with i = 16.0 and 
5'2o = 17.1 mJy. The third row shows the brightest 
radio source, a K4 dwarf with i = 16.7 and S20 = 203 
mJy. The bottom row shows an Ml dwarf with weak 
Ha emission (not visible on the scale of this figure), 
i = 17.8, and 6*20 = 6.6 mJy. 



ing only those matches which would qualify for SDSS 
spectroscopy reduced the random samples to 225 ± 22. 
Applying the su ccess rate of spectral sampling of SDSS 
targets (92.5%; Blanton et al. 20031 ) results in an ex- 
pected random sample size of 208 ± 20. 

We can make an educated guess that the remain- 
ing random samples consist almost entirely of stars: at 
these very bright magnitudes, the SDSS star — galaxy 
separation mechanism is quite effective at differentiat- 
ing between point sources and extended sources. Be- 
sides stars, the most common type of optical point 
source is a quasar, but these are rare at i < 19.1. 
For example, the highest fraction of quasars in SDSS 
can be found at the North Galactic Pole, where there 
are ab out 100 sta,rs for each quasar (at magnitudes 
i < 19: IJuric et al.ll2008l ). Close to the Galactic plane, 
the quasar fraction is much smaller. 

The final step in the selection process was to vi- 
sually classify each source according to its radio mor- 
phology. Only those objects which appeared by eye 
to be unresolved in their FIRST image were retained. 
Selecting these objects results in a random sample size 
of 117 ± 14. Applying the spectral observation success 
rate reduces that value to 110 ± 14. Estimating one 
per cent of the sample are quasars as described above, 
the final estimate is 108 ± 13. This number is essen- 
tially identical to the size of the candidate radio stars 
sample, which consists of 104 stars from DR6 and 8 ad- 
ditional stars with spectroscopy performed later than 
that of DR6. 

The above comparison shows that most or all of the 
potential radio stars are actually chance alignments of 
SDSS stars with unrelated FIRST sources. The vari- 
ance in the random samples is large enough that there 
may be several real radio stars in the candidates sam- 
ple or none at all. This result indicates that radio 
stars are extremely rare or non-existant in the range 
15 < i < 19.1, 5*20 < 1.25 mJy. 

3.7. The Fraction of Radio Stars in the SDSS 

We can use the relative sizes of the candidate and 
random samples, along with the sample completeness, 
to calculate an upper limit on the fraction of radio 
stars in the SDSS. The candidate sample contains 104 
stars with spectra in DR6; 98 of those are in the magni- 
tude range 15 < i < 19.1. As discussed previously, the 
estimated number of contaminating sources is 108±13. 
We can therefore state with 97.5% confidence that 
there are no more than 16 radio stars in the sample 
of candidates (using the one-sided 2-a error estimate) . 
The completeness estimate has the following contribu- 
tions: 1) the completeness of the FIRST surv ey at 1.25 
mJy is approximately 85-90% (jBecker et al.i ,1995): 
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2) the completeness of the radio — optical matching 
within 1" is 96%; 3) the SDSS spectroscopic target- 
ing algorithm selects all 15 < i < 19.1 objects within 
2" of a FIRST source (100% completeness); 4) given 
fiber collisions, the success rate of spectroscopi c ob- 
servations is about 92. 59^ (|Blanton et al.ll2003l ). We 
therefore estimate that the sample is about 75% com- 
plete, and thus that there are no more than 21 radio 
stars in the SDSS DR6 with 15 < i < 19.1. There 
are approximately 18 million SDSS stars in this mag- 
nitude range, which implies that no more than 1.2 out 
of every million stars in the range 15 < i < 19.1 have 
a radio flux of 5*20 > 1.25 mJy. For stars with i ^ 15, 
that corresponds to an upper limit on radio to optical 
flux ratio of 0.34; for the i « 19.1 stars at the faint 
end, the upper limit on radio to optical flux ratio is 
15. 

4. OPTICAL PROPERTIES 

We showed in the previous section that the sample 
is highly contaminated by interloping AGN. However, 
statistical comparisons of the sample with typical stars 
can highlight the most likely actual radio stars. In 
this section, we examine optical properties of the sam- 
ple: photometric colors, distance from the stellar lo- 
cus, spectral type, and magnetic activity. We also use 
an SDSS control sample selected from a strip of sky 1° 
wide in right ascension. The control sample contains 
point sources with 15 < i < 19.1 and r < 20.5 from 
the region 236 < R.A. < 237, -2.5 < dec. < 65; it 
contains just over 160,000 sources. 

4.1. Photometric colors 

Figure |4] presents the distribution of the candi- 
date radio stars in optical color — color space, com- 
pared to the SPSS stellar locus as parameterized by 



Covey et al.l (|2007l ). The majority of the candidate 
stars lie on the stellar locus. Several stars, how- 
ever, appear to lie along the white-d warf — M dwarf 
(WD-l-dM) bridge (jSmolcic et al.ll2004l ): u-g <2,g~ 
r > 0.3, r — i > 0.7. Such close binary pairs are found 
to be more active than their si n gle field counterparts 
fe.g-. ISilvestri et al.ll2005ll2006l ). ISmolcic et all (|2004[) 
limited their analysis to stars with u < 20.5 to elimi- 
nate those with poor photometry. In the SDSS control 
sample, less than 0.1% of stars lie on the WD-(-dM 
bridge. The candidate radio stars sample contains 
seven stars which lie along the WD-fdM bridge, corre- 
sponding to a much higher fraction. We note however 



^'^The success rate is not biased with respect to the radio sample; 
we verified that approximately 92.5% of our selected sources 
which should have passed the targeting selection do in fact have 
SDSS spectra. 



that all of our WD-|-dM candidates have u-band mag- 
nitudes > 20.5; therefore u-band photometry may have 
artificially moved these sources blueward of the M star 
locus (upper right corner of the stellar locus). None 
of their spectra suggest the presence of a white dwarf 
companion. We note that the star with u — g < is not 
on the WD-|-dM bridge, as it has r — i ^ 0.15. How- 
ever, the spectrum and photometric image do suggest 
that it is a (physical or optical) binary system: a K3 
star with an added blue component. 

A quantitative method of finding outliers from the 
stellar locus , taking phot o metri c errors into account, is 
outlined by ICovey et al.l (|2007l ) . They parameterized 
the stellar locus by finding its 1-a width in the stan- 
dard SDSS colors {u — g , g — r,r — i,i — z) an d Two- 
Micron All Sky Survey ( Skrutskie et al. 2006f) colors 
as a function oi g — i. They chose g — i as the fiducial 
color because it samples the largest wavelength range 
possible without relying on the shallower u ov z mea- 
surements. Using only the SDSS colors, we define a 
"four-dimensional color distance" (4DCD; analogous 
to the seven-dimensional version discussed in CO'iiLj), 
which describes the statistical significance of the dis- 
tance in color space between a target object and the 
point on the stellar locus with the same g — i color as 
the target. The 4DCD is defined by 



4DCD 



E 

.fc=0 



targ 
k 



locus 



a\ (locus) -I- (T^ 



1/2 



(1) 



where Ag — u — g, Xi — g — r, etc. The width of the 
stellar locus is ax (locus) , which refers to the FWHM 
of the locus in color Xk at the same g — i as the target 
object. The error in the target's color, ax, is calculated 
by adding in quadrature its photometric errors in the 
two appropriate filters. 

We characterize outliers from the stellar locus as 
those with a value of 4DCD greater than two, shown 
as circles in Figure HI Extreme outliers, with 4DCD> 
3, are shown as triangles. Three of the stars with 
WD-|-dM colors have 4DCD> 2; these are strong can- 
didates for further investigation. The lower right panel 
compares the cumulative distributions for the candi- 
date radio stars {full line) and the SDSS control sam- 
ple; the two distributions are similar. The largest 
4DCD value of ~ 14 belongs to the possibly binary 
star mentioned earlier, having u — g < Q. The other 
most extreme outlier has 4DCD ~ 6.4. The spectrum 
of that object also shows some evidence that the source 
is part of a multiple system, such as the super-position 



^^Our definition of diff'ers from C07 by a square root operator, 
such that our value is in the units of a color distance. 
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of a late-type with an early-type M star. All other val- 
ues of 4DCD are < 4.3. 



4.2. Spectral type and activity fraction 



Because non-thermal radio emission is a signal of 
activity, we expect that radio emission may correlate 
with strong spectral lines, such as Ha, which are also 
known to signal magnetic activity. We investigate the 
fraction of active stars in our sample; if the sample 
contains some real radio stars which are active, we 
may see an increase over the active fraction of all stars 
(without selecting for radio emission). Previous stud- 
ies have found that the fraction of active M dwarfs is 



a str ong function of spectral type (jWest et al.l 12004 . 
20081 ). tending to increase toward later subtypes with 
a peak around M8 dwarfs. 

An Ha equivalent width was measured for each stel- 
lar S£ectrum_usingthe Hammer f ij3.3p . As discussed 
by IWest et all (|2008l ). the accuracy of such measure- 
ments has been tested via Monte Carlo simulations to 
ascertain how well line strength can be determined at 
a given S/N level. The Ha emission can be recov- 
ered over 96% of the time for all spectral types. Fig- 
ure [5] shows the spectral types of the candidate radio 
stars sample {solid line), compared with all SDSS DR6 
stars with spectra (just under 1 million stars; dotted 
line). Spectral types for the SDSS sample come from 
the automated version of the Hammer, while spec- 
tral types for the candidate radio star sample were 
visually-confirmed ('i )3.3p . Automated Hammer clas- 
sifications are typically accurate to within ±4 sub- 
types for A- G stars, and withi n ±2 subtypes for K 
and M stars ( Covev et al.ll2007h . All spectra for the 
radio stars sample have S/N ratio greater than 2.9. 
As shown in Figure O all of the active candidate ra- 
dio stars are M dwarfs. This is not a s urprising re- 



sult given that most stars a re M dwarfs (jCovev et al, 
2008 : Bochanski et al. 2009f l. and the majority of ac- 
tivity in main seq ue nce stars is seen in M dwarfs (e.g. 



Gizis et al.ll2002l) . IWest et al.l (|2008f ) discussed the 



fraction of active M dwarfs in SDSS DR5 as a function 
of spectral subtype, after removing WD-fdM bridge 
stars from their sample. They found that M0-M3 stars 
have an active fraction of roughly 5-20%, and that the 
fraction increases strongly for M dwarfs of later sub- 
type. Our results for the candidate radio stars are 
consistent with the results of that study for all stars; 
we do not see a significantly higher fraction of active 
stars in the sample of candidate stellar radio sources. 

Figure [5] shows color — color diagrams of the candi- 
date radio stars and the SDSS control sample, with 
the six active stars indicated by triangles. Three of 
the seven stars with WD-|-dM bridge colors ( §4.ip are 
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Fig. 4. — Outliers from the stellar locus. Panels A- 
C: Color — color diagrams with potential radio stars 
plotted as symbols. The red plus signs indicate can- 
didate radio stars consistent with being on the stellar 
locus, having 4DCD < 2 (see text and Eq. 1 for defini- 
tion); cyan circles indicate those with 2 < 4DCD < 3; 
green triangles indicate those with 3 < 4DCD. The 
dashed-line rectangle in panel A surrounds the stars 
with WD-hdM bridge colors. The SDSS stellar lo- 
cus is shown by t he solid lines, as parametrized by 



Covey et al.l ([20071). The three lines indicate the po- 
sition of the stellar locus and its interquartile width, 
projected onto two dimensions. Panel D: Cumulative 
distributions of 4DCD for the candidate radio stars 
sample (solid line) and for the SDSS control sample 
(dotted line). 
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active; two of those are outliers from the steUar lo- 
cus (4DCD > 2). T his result is in agreement with 
Silvestri et aD (|2006h . who suggest that 20-60% of aU 
WD+dM binaries are magnetically active. The lower 
right panel of Fig. [5] shows the g~i distribution of the 
SDSS control sample {dotted line), the candidate radio 
stars (thin solid line), an d the active stars (thick solid 
line). ICovey et al. ( 20071 ) showed that g — i correlates 
strongly with stellar spectral type. The distribution 
is bimodal in flux limited samples: it is biased toward 
red stars, which are the most common, but also blue 
stars, which can b e see n from much grea. t er dist ances- 
West et alJ (|2004 ) and iBochanski et all (|2007t ) com- 
pared the colors of active to inactive M dwarfs in the 
field (i.e., not part of a binary system). They found 
no significant differences between the two populations. 
The current radio stars study includes three active field 
M dwarfs (i.e., M dwarfs which do not have WD-fdM 
bridge colors). Although the sample is too small to 
make a strong sta t emen t in c omparison with the stud - 
ies of I West et all (|2004l) and IBochanski et all (|2007h . 
we note that our results for stars with radio emission 
are consistent with their conclusions for all stars. 

5. RADIO PROPERTIES 
5.1. Radio Spectral Slope 

Using the multiple-wavelength radio catalog pre- 
sented by KI08, we find that eight of the candidate 
radio stars were detected in a second radio sky sur- 
vey: either at 92 cm in WENSS or at 6 cm in GB6 (or 
both). We can therefore determine the radio spectral 
index a (where (X v"), and we report those results 
here. We compare the results with samples of the type 
of AGN that may be contaminating the sample ( §3.6p . 

The value of the spectral index is a clue about the 
environment at the source of emission. Non-thermal 
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Fig. 6. — Activity in the candidate radio stars sam- 
ple. Panels A-C: Color — color diagrams with the in- 
active stars indicated by red plus signs and the active 
stars (showing Ha in emission) indicated by blue trian- 
gles. The dashed-line square in panel A surrounds the 
stars with WD-hdM bridge colors. The SDSS stellar 
locus is shown by the solid lines, as parametrized by 
Covey et al.l (|2007l) . The three lines indicate the po- 
sition of the stellar locus and its interquartile width, 
projected onto two dimensions. Panel D: The thin 
solid line shows the distribution oi g — i for the can- 
didate radio stars with i < 19.1. The thick line indi- 
cates the six active stars. The dotted line corresponds 
to stars in the SDSS control sample with photometric 
errors less than 0.1 mag in g and i. 
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Fig. 5. — Spectral types of the candidate radio stars 
(solid line). The thick line indicates the six stars with 
Ha in emission. SDSS stars with 15 < i < 19.1 and 
r < 20.5 are also shown (dotted line). 
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radio emission is typically due to synchrotron (rel- 
ativistic) or gyrosynchrotron (semi-relativistic) elec- 
trons accelerating in a magnetic field. Synchrotron and 
gyrosynchrotron processes result in a negative spectral 
slope (a ~ —0.8) in the optically-thin case, and a flat 
or positive slope in the optically-thick case. Such emis- 
sion has been detected from M dwarfs in both their 
flaring (rapi dly-variable) and qui e scent (non-fiaring 



states (e.g.. iGiidel fc Benz 



Bastian fc Bookbinder 11987 



1996t iLarge et al.l \l98 



Osten et al 



20061) 



Fifty stars in our sample lie within the WENSS sky 
coverage and 106 lie within the GB6 sky coverage. 
Following KI08, we use 30" as the FIRST— WENSS 
matching radius and 70" as the FIRST — GB6 match- 
ing radius. These choices result in estimates of 99% 
completeness with 92% efficiency (WENSS), and 98% 
completeness with 79% efficiency (GB6; see Table 2 of 
KI08). Six sources have 92 cm detections and three 
have 6 cm detections; only one source was detected 
at all three wavelengths. FIRST is much deeper than 
WENSS or GB6; therefore only sources which are very 
bright or have very steep spectral slopes can be de- 
tected in more than one survey. Despite the poorer po- 
sitional accuracy of GB6 and WENSS, these matches 
are highly reliable because of the their lower source 
sky densitvF^ Figure [7] presents the radio spectra of 
these sources. Four of the stars are M dwarfs, two of 
which are active ( §4.2p . 

From the KI08 catalog, we selected a sample of 
AGN using the radio criteria that were applied to the 
radio stars sample: 5*20 < 1.25 mJy, unresolved in 
FIRST (using Eq. 3 of KI08), and a WENSS match 
within 30". These criteria select the type of objects 
that could be contaminating the sample of potential 
radio stars discussed in this section. The stellar sam- 
ple has a similar spectral slope distribution to the AGN 
sample, whose median (mean) spectral index is -0.67 
(-0.54). The spectral slope distributions do not allow 
a definitive conclusion as to whether any of the can- 
didate radio stars with multiple radio detections are 
actual radio stars or are contaminating AGN. 

5.2. Variability 

It is possible that the radio fluxes varied over the 
course of observations by the different surveys. Both 
stellar radio sources and AGN are known to vary in the 
radio. For example. iBerge 2 (|2002h observed variability 
in the non-flaring radio emission from cool dwarfs over 
timescales of merely hours. (However, the fluxes in 
that study were far below the flux limit of the current 



study.) The likely AGN contaminants are quasar core 
sources, which are known to v ary with timescales rang- 
ing from days to ye ars (e.g., Rvs fc Machalskil 1990l : 
Barvainis et al.ll2005l ) . About half of the candidate ra- 
dio stars sample was observed muliple times in FIRST, 
most with timescales of less than a week. None showed 
significant variability at 20 cm. We note that nearly 
half of the radio stars found by II99 varied in FIRST 
at the 4ct level. 

6. CONCLUSIONS 

We performed a search for radio stars by combining 
radio and optical observations from FIRST (20 cm) 
and SDSS. This is the first large-scale search for ra- 
dio stars using an optical survey faint enough to in- 
clude a large number of M dwarfs. Many of these 
late-type stars are known to be magnetically active. 
A sample of 112 candidates was selecting using the 
following criteria: optical point source morphology as 
determined by the SDSS photometric pipeline, radio 
and optical positions matched within 1", an optical 
magnitude r < 20.5 to ensure reliable determination 
of optical morphology, radio flux ^20 > 1-25 to elimi- 
nate spurious sources, a spectrum visually classified as 
stellar, and radio point source morphology. We esti- 
mated sample contamination using random matching 
and similar selection criteria. The size of the random 
samples (108 ± 13) suggests that the potential radio 
stars are heavily contaminated by optically-faint ra- 
dio quasars in chance alignment with a foreground 
star. The main ambiguity in determining radio — 
optical matches stems from uncertainty as to whether 
the star is actually the source of the radio emission. It 
may be possible to overcome this problem with care- 
ful follow-up using very long baseline interferometry 
to determine proper motions of the radio sources. 

In we calculated the upper limit on the fraction 
of radio stars (with 820 > 1.25 mJy) no more than 1.2 
per million stars in the magnitude range 15 < i < 19.1. 
We note that some M stars have been observed to oc- 
casionally flare brightly in the radio; however, these 
stars have a very small duty cycle and very few will 
be detected in a single-epoch survey such as FIRST. 
While some M dwarfs have been shown to have con- 
stant radio emission (e.g., Bergei 20021 : Osten et al 



^^For each object, we verified that the radio star candidate is the 
nearest FIRST source to the GB6 or WENSS match. 



20061 ) ■ it is at much fainter levels than radio emission 
in their flare state by about an order of magnitude. 
Our results effectively rule out a population of radio 
bright late-type stars. 

We compared the candidate sample to SDSS stars 
in the same magnitude range and investigated their 
distributions in optical color, stellar type, magnetic 
activity, and distance from the stellar locus in color — 
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Fig. 7. — Le/t; Radio spectra of the candidate radio 
stars detected at more than one radio wavelength. Di- 
amonds indicate inactive stars, asterisks indicate ac- 
tive stars. The visuaUy-confirmed spectral type and 
the power law index are labeled for each source. Right: 
distribution of FIRST fluxes for the 44 candidate ra- 
dio stars which lie in the region of the sky observed by 
WENSS, but were not detected at 92 cm. 
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Fig. 8. — Parameter space probed by combinations of 
radio and optical surveys. Stellar radio sources de- 
tected by H99 are indicated by symbols. Dashed lines 
show constant radio-to-optical flux ratios of 1:1000, 
1:100, 1:1, 100:1 (from left to right). The FIRST— 
SDSS matching corresponds to the current paper. 
LSST and EMU are future southern sky surveys. 



color space. The two data sets show similar distri- 
butions. However, there is a higher fraction of stars 
with WD+dM bridge colors among the candidate sam- 
ple. Three of the radio star candidates with WD-|-dM 
bridge colors are significantly offset from the stellar 
locus in the four-dimensional color space (4DCD> 2). 
Two of those three stars are magnetically active; a 
third star on the WD-f dM bridge is also active. These 
four stars are good candidates for continuing investi- 
gation. 

We searched for radio detections at two other wave- 
lengths, 6 cm and 92 cm. Eight sources were detected 
at more than one radio wavelength. The spectral index 
distribution for the stellar sources is similar to the dis- 
tribution for a sample of possible AGN contaminants. 
Therefore the subset of stellar radio candidates may 
be plagued by the same AGN contamination seen in 
the overall sample. 

This study shows that FIRST and SDSS are not 
a good pair of surveys for the study (or discovery) 
of radio stars: stars bright enough at 20 cm to ap- 
pear in FIRST are probably above the m = 15 satu- 
ration limit of the SDSS. Figure [8] shows the radio — 
optical parameter space probed by the II99 study and 
the FIRST — SDSS correlation presented in this pa- 
per. FIRST is currently the largest deep radio survey 
available; the figure indicates that much fainter radio 
data is required in order to discover more stars of the 
population found by H99. A much deeper sky sur- 
vey at the same frequency as FIRST is likely to be 
carried out in the southern hemisphere using the Aus- 
tralian Square Kilom eter Array PathfindeiLj (ASKAP; 



Johnston et"al]l2009f ). One of ASKAP 's invited pro- 



posals involves a project known as "EMU: Evolution- 
ary Map of the Universe" (R. Norris et al. in prepa- 
ration), which includes a large southern sky survey 
down to 50 /iJy {barms) at 20 cm. Although the 
main science driver for the ASKAP EMU project is 
the study of AGN and galaxy evolution, the field of 
radio stars will benefit immensely by cross-correlating 
the EMU survey with an all-sk y optical survey su ch as 
the Guide Star Catalog fGSC: llasker et al.lll99dl ) or a 
large southern sky survey s uch as the Large S ynoptic 
Survey Telescop' 



ky survey su 
S (LSST;|i 



vezic et al.ir2008bh . 



The GSC, the faintest of the optical surveys used 
in the H99 study, covered the sky in the approximate 
magnitude range of 7-15. As shown in Fig. [51 the com- 
bination of the GSC and EMU will be sensitive to stars 
with radio-to-optical flux ratios from less than 1:1000 
to those with the highest flux ratios found by H99. 
Covering half the sky and extending 20 times fainter 



^•^ http://www. atnf. csiro.au/projects/askap/ 
-'^ http://www. Isst. org 
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in the radio than FIRST, GSC— EMU should result in 
the largest sample of candidate radio stars to date. 

The LSST is a multi-epoch optical survey which 
will observe a quarter of the sky every three nights 
and detect point sources down to r ~ 27 (r ~ 24.5 
for a single exposure). The radio — optical parameter 
space covered by a potential LSST — EMU matching 
is also shown in Fig. \8\ The H99 survey found only a 
couple of radio stars of the type to which LSST — EMU 
will be sensitive. The EMU survey, 20 times fainter 
than FIRST, will probe a 20-90 times larger volume 
in the Milky Way disk, and thus may find on the order 
of 100 new sources of faint stellar radio emission. In 
addition, the LSST will be able to recognize active M 
dwarfs via their UV flaring. This method could be 
more efficient than looking for Ha emission in SDSS 
spectra, which are available for only about 0.25% of M 
dwarfs detected by the SDSS. The advent of the EMU 
survey to radio astronomy, combined with large optical 
surveys, is likely to increase the number of known radio 
stars by orders of magnitude. 
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Table 1 

Properties of the radio star candidates sample 





FIRST'-' 


separation^ 


distance^ 






5pcak 


S20 




SDSS model magnitudes 






catalog 


r,.a. 


dec. 


r,.a. 


dec. 


["] 


[pc] 


typc*^ 


active?'^ 


[mJy /beam] 


[mJy] 


u 


g 


r 


i 


z 


4DCD'' 


ID*= 


9.7469648 


0.31168685 


9.74689 


0.31172 


0.29 


319.640 


M2 


no 


1.93 


2.01 


21.3 


18.8 


17.4 


16.4 


15.9 


1.7 


55105 


27.6658730 


-1.13948400 


27.66569 


-1.13928 


0.99 


397.880 


M2 


no 


1.82 


1.96 


21.6 


18.9 


17.5 


16.6 


16.2 


1.8 


155708 


44.5010990 


1.22593510 


44.50115 


1.22574 


0.73 


9830.540 


F8 


no 


11.90 


12.40 


19.2 


18.3 


18.1 


18.0 


18.0 


2.8 


253325 


112.0061000 


38.06249000 


112.00596 


38.06232 


0.73 


598.000 


Ml 


no 


1.44 


1.60 


22.0 


19.8 


18.3 


17.5 


17.1 


1.9 


612859 


114.5806600 


48.39231500 


114.58059 


48.39211 


0.76 


2750.950 


K2 


no 


1.27 


1.81 


20.3 


18.7 


18.1 


17.9 


17.8 


0.8 


634578 


117.5012900 


34.98278400 


117.50105 


34.98291 


0.85 


691.950 


Ml 


no 


1.83 


1.56 


21.3 


18.9 


17.6 


17.0 


16.6 


2.5 


660499 


117.5527500 


24.83105600 


117.55283 


24.83080 


0.96 


1139.890 


MO 


no 


2.97 


2.85 


23.1 


19.7 


18.4 


17.9 


17.6 


2.2 


661031 


118.5574000 


39.62222400 


118.55721 


39.62228 


0.58 


2268.550 


F8 


no 


1.27 


1.29 


16.9 


15.9 


15.6 


15.5 


15.5 


0.4 


670267 


118.6814000 


33.29291800 


118.68152 


33.29275 


0.70 


515.780 


Ml 


no 


1.80 


1.93 


21.3 


18.8 


17.4 


16.8 


16.4 


1.1 


671427 


119.3160100 


30.17021500 


119.31614 


30.17024 


0.41 


936.480 


K4 


no 


3.57 


3.82 


19.8 


17.7 


16.7 


16.4 


16.2 


3.5 


677301 


119.7809500 


17.07585700 


119.78081 


17.07571 


0.71 


1613.270 


MO 


no 


2.37 


1.79 


23.2 


20.2 


18.9 


18.4 


18.1 


1.2 


681437 


121.2057400 


53.11331800 


121.20609 


53.11341 


0.83 


1070.120 


M3 


no 


3.81 


4.97 


22.6 


21.1 


19.8 


18.8 


18.3 


3.7 


694900 


121.3575900 


33.92661500 


121.35767 


33.92649 


0.51 


1040.680 


M3 


no 


1.74 


1.25 


23.6 


21.6 


20.2 


19.2 


18.6 


1.2 


696337 


122.0831300 


40.20902400 


122.08297 


40.20927 


0.99 


12403.900 


KO 


no 


1.48 


1.36 


20.2 


19.2 


18.9 


18.8 


18.8 


1.2 


703419 


122.4434000 


15.34121000 


122.44323 


15.34111 


0.69 


795.440 


M2 


yc 


1.45 


3.43 


22.8 


20.2 


18.9 


18.0 


17.5 


1.9 


706934 


122.5364400 


39.29135700 


122.53619 


39.29131 


0.72 


887.850 


K4 


no 


4.99 


4.27 


18.4 


16.6 


15.9 


15.7 


15.6 


1.1 


707887 


123.5006500 


29.96177800 


123.50080 


29.96161 


0.77 


535.690 


M4 


m 


6.21 


5.72 


23.0 


21.6 


20.2 


18.8 


18.1 


2.0 


717543 


123.5311000 


7.57838460 


123.53087 


7.57825 


0.94 


395.760 


M4 


no 


1.01 


1.58 


23.6 


20.8 


19.3 


18.1 


17.4 


2.4 


717871 


123.8783600 


27.77808900 


123.87866 


27.77810 


0.97 


4003.250 


G2 


no 


2.84 


2.45 


19.5 


18.3 


17.9 


17.8 


17.7 


1.9 


721293 


124.0573700 


17.92303000 


124.05750 


17.92286 


0.76 


1131.870 


K4 


no 


1.38 


1.88 


19.2 


17.4 


16.7 


16.4 


16.3 


1.2 


723055 


125.2091100 


42.31718800 


125.20885 


42.31716 


0.69 


192.380 


Ml 


no 


4.34 


4.27 


19.7 


17.1 


15.7 


14.9 


14.5 


0.8 


734833 


126.3005300 


17.33529000 


126.30050 


17.33502 


0.98 


1258.650 


Ml 


no 


2.15 


1.86 


23.1 


20.8 


19.4 


18.7 


18.4 


1.3 


746124 


127.3762300 


47.77290800 


127.37631 


47.77266 


0.92 


1149.710 


K4 


no 


194.00 


203.00 


19.9 


17.9 


17.1 


16.7 


16.6 


1.4 


756772 


128.4548800 


28.86214000 


128.45489 


28.86231 


0.61 


315.370 


M3 


no 


9.40 


9.77 


22.1 


19.7 


18.4 


17.1 


16.5 


2.4 


767820 


129.9784600 


17.21355300 


129.97861 


17.21342 


0.69 


483.870 


M4 


yc 


2.50 


2.46 


23.2 


21.6 


20.2 


18.8 


18.0 


2.4 


783928 


131.7995200 


8.87368060 


131.79954 


8.87363 


0.20 


1409.130 


K3 


no 


1.67 


1.89 


19.9 


17.9 


17.2 


16.9 


16.8 


1.2 


803383 


132.3854500 


39.56273500 


132.38580 


39.56271 


0.98 


638.660 


MO 


no 


1.43 


2.29 


21.6 


18.9 


17.6 


17.0 


16.7 


1.0 


809871 


133.7192300 


22.38336800 


133.71919 


22.38349 


0.46 


3791.930 


F8 


no 


3.15 


3.00 


18.2 


17.2 


16.9 


16.8 


16.8 


0.3 


824171 


133.9574700 


16.77287300 


133.95775 


16.77288 


0.95 


2133.210 


G3 


no 


16.20 


17.10 


17.6 


16.5 


16.1 


16.0 


16.0 


2.9 


826748 


138.3544300 


16.92733700 


138.35438 


16.92730 


0.22 


930.940 


M2 


no 


4.84 


4.56 


22.5 


20.4 


19.1 


18.3 


17.8 


2.6 


874572 


139.0118400 


30.04233700 


139.01158 


30.04221 


0.93 


541.220 


M4 


ye 


2.16 


2.39 


23.2 


21.4 


19.9 


18.7 


18.0 


1.3 


881748 


139.3560700 


10.26608300 


139.35585 


10.26599 


0.85 


1857.790 


MO 


no 


5.43 


5.64 


23.1 


20.7 


19.4 


18.9 


18.5 


1.6 


885740 



Note. — The table of 112 candidate radio stars is available at http://www.astro.washington.edu/users/akirnball/radiocat/radiostars/. 
^Right ascension and declination are given in decimal degrees. 
^Offset between FIRST and SDSS positions. 

'^Distance was determined using the photometric parallax relation of llvezic et al.| ( |2008aD ; the relation is valid only for stars on the main sequence. 
■^Visually-confirmed spectral classification. 

"yes" indicates a spectrum with reliable Hq emission. 
^Four-dimensional color distance from the stellar locus, defined in §4.1 of the text. 
^Internal ID of the source in the radio catalog of KI08. 



Table 2 

Candidate radio stars" with resolved radio emission 





FIRST'' 


separation^ 


distance^ 






Spoak 


S20 




SDSS model magnitudes 




catalog 




r,.a. 


dec. 


r,.a. 


dec. 


["] 


[pc] 


typc*^ 


active?^ 


[m Jy /beam] 


[mJy] 


u 


g 


r 


i 


z 


4DCD^ 


ID« 


123.1976100 


20.64439200 


123.19768 


20.64418 


0.80 


351.420 


Ml 


no 


1.25 


2.11 


20.6 


18.0 


16.6 


15.9 


15.5 


0.6 


714327 


129.6339900 


13.72514400 


129.63398 


13.72527 


0.45 


2147.920 


G4 


no 


16.60 


30.00 


18.4 


17.2 


16.7 


16.6 


16.5 


4.6 


780377 


159.1123200 


12.80009800 


159.11230 


12.80003 


0.26 


12229.900 


F7 


no 


68.30 


103.00 


20.2 


19.3 


19.1 


19.0 


18.9 


3.5 


1101138 


179.5614700 


37.07555200 


179.56123 


37.07574 


0.97 


6233.260 


K2 


no 


1.39 


2.07 


21.7 


20.4 


19.8 


19.6 


19.4 


1.6 


1324909 


179.9687500 


44.91827900 


179.96864 


44.91834 


0.36 


6285.120 


G5 


no 


3.73 


6.41 


20.4 


19.5 


19.1 


18.9 


18.8 


3.1 


1329461 


183.0861100 


27.54304500 


183.08620 


27.54281 


0.89 


687.640 


M2 


no 


4.66 


7.81 


23.2 


20.3 


18.8 


18.0 


17.5 


0.8 


1363672 


189.3334100 


53.88457500 


189.33319 


53.88464 


0.52 


2723.890 


Gl 


no 


3.37 


6.85 


18.8 


17.7 


17.2 


17.1 


17.0 


2.2 


1431967 


198.7216200 


10.62265900 


198.72157 


10.62287 


0.78 


1404.190 


M2 


no 


3.75 


7.46 


22.8 


21.2 


19.9 


19.1 


18.4 


5.6 


1534251 


217.8390300 


8.53663970 


217.83908 


8.53640 


0.88 


2057.150 


K4 


no 


1.92 


4.61 


21.1 


19.1 


18.3 


17.9 


17.8 


1.0 


1742370 


228.9030100 


6.33068160 


228.90314 


6.33059 


0.56 


2877.940 


Gl 


no 


2.90 


5.90 


18.3 


17.4 


17.0 


16.9 


16.8 


3.0 


1863654 


240.3610400 


9.07629920 


240.36122 


9.07639 


0.72 


1297.210 


Ml 


no 


30.30 


62.30 


23.4 


20.5 


19.2 


18.5 


18.2 


0.5 


1985947 


241.6174600 


32.45298600 


241.61725 


32.45302 


0.65 


934.290 


Ml 


no 


13.60 


26.00 


22.1 


19.7 


18.5 


17.7 


17.3 


3.4 


1998928 


251.0205900 


26.75291400 


251.02068 


26.75277 


0.60 


9124.150 


F8 


no 


3.67 


7.14 


19.9 


18.9 


18.7 


18.5 


18.5 


2.8 


2089662 


253.7596800 


32.12100800 


253.75963 


32.12106 


0.24 


1440.760 


K4 


no 


1.02 


4.60 


21.2 


19.1 


18.0 


17.6 


17.4 


1.7 


2113922 


254.7470100 


23.40834500 


254.74695 


23.40835 


0.21 


548.530 


MO 


yc 


1.27 


2.70 


20.7 


18.6 


17.3 


16.6 


16.2 


3.9 


2122609 


258.2688900 


32.22357000 


258.26875 


32.22338 


0.81 


2846.860 


F9 


no 


1.44 


3.51 


18.2 


17.2 


16.9 


16.7 


16.7 


3.0 


2152684 



Note. — The tabic of 16 sources with resolved radio emission (which passed all other selection criteria) is available at http://www.astro.washington.edu/users/akim.ball/radiocat/radiostars/. 
'^Right ascension and declination are given in decimal degrees. 
^Offset between FIRST and SDSS positions. 

^Distance was determined using the photometric parallax relation of llvezic et al.| ||2008aD ; the relation is valid only for stars on the main sequence. 
'^Visually-confirmed spectral classification. 

"yes" indicates a spectrum with reliable Hq emission. 
^Four-dimensional color distance from the stellar locus, defined in §4.1 of the text. 
^Internal ID of the source in the radio catalog of KI08. 



Table 3 

"Candidate radio stars" with complex radio emission 



SDSS'^ 




separation^ 


distance'^ 






•Speak 


S'20 




SDSS model maj 


^nitudes 




catalog 




R. A. 


dec. 


R. A. 


dec. 


["] 


[pc] 


typc^ 


active?^ 


[mjy/beam] 


[mJy] 


u 


g 


r 


i 


z 


4DCD' 


ID" 


8.8779417 


-10.33142200 


8.87798 


-10.33148 


0.25 


365.070 


M3 


no 


4.87 


5.79 


22.6 


19.7 


18.3 


17.2 


16.6 


0.7 


50158 


16.4108070 


0.04541176 


16.41080 


0.04553 


0.43 


1605.120 


K7 


no 


1.36 


4.61 


22.9 


20.3 


19.0 


18.5 


18.3 


3.2 


92777 


30.9980570 


-9.00080110 


30.99826 


-9.00074 


0.75 


1409.850 


G2 


no 


1.41 


1.97 


16.7 


15.6 


15.2 


15.1 


15.1 


0.4 


174929 


31.6274380 


-8.36097310 


31.62769 


-8.36088 


0.96 


7284.830 


F8 


no 


1.05 


1.96 


20.0 


19.0 


18.6 


18.5 


18.5 


1.3 


178794 


115.4349200 


33.59706800 


115.43505 


33.59718 


0.57 


660.370 


M3 


no 


117.00 


137.00 


23.3 


20.8 


19.3 


18.3 


17.8 


0.5 


642007 


116.5239800 


33.14653100 


116.52415 


33.14652 


0.51 


1299.940 


G5 


no 


24.40 


42.50 


17.2 


16.0 


15.6 


15.5 


15.4 


0.8 


651724 


118.6515700 


31.04799200 


118.65183 


31.04788 


0.89 


976.060 


K2 


no 


3.14 


5.70 


18.2 


16.6 


15.9 


15.7 


15.6 


1.1 


671131 


120.2777900 


31.98149700 


120.27754 


31.98134 


0.95 


3123.470 


G5 


no 


1.40 


1.80 


18.7 


17.5 


17.2 


17.0 


17.0 


1.1 


686143 


122.4168500 


9.98404330 


122.41704 


9.98402 


0.69 


2606.500 


G8 


no 


9.55 


22.70 


19.4 


18.1 


17.6 


17.4 


17.3 


1.1 


706655 


125.7096000 


12.55254200 


125.70952 


12.55259 


0.33 


857.030 


M3 


ye 


7.58 


15.40 


23.0 


21.4 


19.9 


19.0 


18.5 


4.2 


740037 


128.9825100 


32.53530500 


128.98271 


32.53522 


0.69 


919.880 


F8 


no 


29.40 


53.70 


23.7 


22.6 


21.0 


19.9 


19.3 


3.5 


773185 


131.5889500 


13.51593600 


131.58898 


13.51595 


0.10 


881.950 


M3 


m 


4.53 


6.21 


23.3 


21.4 


20.0 


18.9 


18.4 


1.3 


801062 


131.9006900 


15.72048400 


131.90071 


15.72021 


0.99 


916.410 


M3 


no 


1.08 


3.20 


22.6 


21.5 


20.0 


19.1 


18.6 


3.7 


804509 


132.1393200 


3.49329870 


132.13913 


3.49330 


0.69 


7584.010 


M2 


no 


5.51 


13.10 


20.3 


19.4 


19.0 


18.9 


18.9 


2.7 


807208 


133.3787500 


3.04040040 


133.37881 


3.04044 


0.27 


1017.230 


Ml 


no 


5.67 


8.47 


22.9 


20.4 


19.1 


18.3 


17.9 


0.7 


820636 


134.9379600 


15.68712000 


134.93808 


15.68694 


0.77 


1733.890 


K4 


no 


29.40 


32.10 


21.0 


19.0 


18.1 


17.7 


17.6 


1.6 


837341 


135.5234500 


-0.30775603 


135.52351 


-0.30753 


0.84 


1238.090 


K7 


no 


2.59 


3.77 


21.8 


19.4 


18.2 


17.7 


17.5 


1.1 


843673 


140.7074600 


0.98319174 


140.70767 


0.98337 


0.98 


2584.960 


K2 


no 


10.50 


16.70 


19.6 


18.3 


17.7 


17.5 


17.4 


1.6 


900871 


141.8337800 


14.43029200 


141.83377 


14.43037 


0.28 


392.550 


M4 


no 


10.60 


11.40 


22.6 


20.4 


19.0 


17.7 


17.1 


1.3 


912857 


147.9741100 


0.86803458 


147.97390 


0.86804 


0.76 


426.410 


M4 


ye 


5.80 


7.25 


22.7 


20.6 


19.2 


17.9 


17.2 


2.1 


979415 


148.9287500 


40.24028700 


148.92868 


40.24017 


0.46 


301.060 


M4 


no 


2.90 


5.75 


22.6 


20.2 


18.8 


17.4 


16.7 


1.2 


989703 


153.7401400 


14.70054700 


153.74037 


14.70052 


0.80 


1375.330 


K7 


no 


7.83 


11.60 


21.7 


19.2 


18.1 


17.6 


17.4 


1.6 


1042751 


154.0058700 


12.36168900 


154.00611 


12.36181 


0.95 


2299.170 


K7 


no 


7.56 


8.66 


22.4 


20.2 


19.1 


18.7 


18.5 


2.0 


1045679 


160.5310600 


3.99054340 


160.53126 


3.99068 


0.88 


809.980 


M3 


no 


5.64 


8.20 


23.0 


21.3 


20.0 


18.9 


18.2 


2.6 


1116725 


172.0590700 


21.93210400 


172.05921 


21.93225 


0.71 


602.870 


M4 


no 


4.60 


6.45 


23.9 


21.6 


20.3 


18.9 


18.2 


2.3 


1242017 


178.5105200 


3.21151700 


178.51063 


3.21162 


0.54 


10917.700 


F4 


no 


54.30 


72.20 


19.9 


19.0 


18.7 


18.7 


18.7 


1.1 


1313182 


186.8444200 


46.53642400 


186.84434 


46.53632 


0.42 


2233.410 


MO 


no 


8.41 


14.90 


22.5 


20.6 


19.5 


18.9 


18.5 


5.0 


1404584 


190.8415000 


46.59573600 


190.84160 


46.59557 


0.64 


2809.610 


G5 


no 


1.06 


3.81 


18.9 


17.8 


17.3 


17.2 


17.1 


2.6 


1448337 



Note. — The table of 60 sources with complex radio emission (which passed all other selection criteria) is available in the electronic version of this paper, and is also downloadable from 
http : / /www . astro. Washington, edu/users/akimball/radiocat/radiostars/. 

^Right ascension and declination arc given in decimal degrees. 

^Offset between FIRST and SDSS positions. 

■^Distance was determined using the photometric parallax relation of llvezic et al.| d2Q08a[) : the relation is valid only for stars on the main sequence. 
'^Visually-confirmed spectral classification. 

"yes" indicates a spectrum with reliable Ha emission. 
^Four-dimensional color distance from the stellar locus, defined in §4.1 of the text. 
^Internal ID of the source in the radio catalog of KI08. 



